GENERAL DECOMPOSITION OF RADIAL 

FUNCTIONS ON M" AND APPLICATIONS TO A^-BODY 

QUANTUM SYSTEMS 

P^ ! CHRISTIAN HAINZL^ AND ROBERT SEIRINGER^ 

Cn ' Abstract. We present a generalization of the FefTerman-de la Llave 

r^ . decomposition of the Coulomb potential to quite arbitrary radial func- 

^ ' tions V on R" going to zero at infinity. This generalized decomposition 

^~~i I can be used to extend previous results on A'^-body quantum systems 

QQ ■ with Coulomb interaction to a more general class of interactions. As an 

^— ( ' example of such an application we derive the high density asymptotics 

, of the ground state energy of jellium with Yukawa interaction in the 

Cn ' tjiermodynamic limit, using a correlation estimate by Graf and Solovej 

> ■ 

o. 

O I 1. Introduction 

^^ ■ For the description of physical systems interacting with Coulomb forces, 

,^ . the Fefferman-de la Llave decomposition of the Coulomb potential has pro- 

Mh| ved very useful. It was introduced in ||FL|| , where it was used in the proof of 

^^ • stability of relativistic matter. It states that, for x S M^, 

c3 ■ 1 1 f°° 1 ^ ^ 

\x\ vr Jn r^ 

• ■ 

^^ ' where Xr{x) = 0{r — \x\) is the characteristic function of a ball of radius 

S^ . r centered at the origin, and * denotes convolution on M^. Except for the 

H \ constant l/vr, it is easily checked that (|l|) holds true, since the right side is 

■ - - ' a radial, homogeneous function of order —1. 

We are interested in a generalization of (|l|) to arbitrary radial functions 
V on M" going to zero at infinity, with a weight function g{r) replacing 1/r^ 
in the integrand (see (|^). In Theorem |l] below we present a simple and 
straightforward derivation of such a decomposition, assuming V to satisfy 
some decrease and regularity properties specified below. We give an explicit 
expression for the weight function g{r), which turns out to be related to 
the [n/2\ + 2'th derivative of V. Here [•] denotes the GauB bracket, i.e., 
[m\ = max{n € No,n < m}. In particular, the case of a positive g is of 
interest, implying positive definiteness of the function V . 
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Decompositions of the type @ have also been studied in [Gl|] for the case 
of positive and bounded functions V, where the decomposition (^) is referred 
to as "scale mixtures of Euclid's Hat" . The improvement of our Theorem |l| 
compared to |G1| lies in its concise proof and the specific, unified form of 



the formula for g. 

The decomposition we derive in the next section is of particular interest 
in generalizing results on A^-body quantum systems so far only applicable 
to the Coulomb potential. For instance, it can be used in estimating the 
"indirect part" of the interaction energy, as was done in the Coulomb case in 



[L0|, or to determine the validity of Hartree-Fock approximations (see 



and |GS|). This is due to the fact that with the help of (y) the interaction 

(2) 



energy in an A^-particle quantum system can be decomposed as 

r /"OO 



where X"^'^ denotes the one-particle multiplication operator Xri/ ~ z)i being 
actually a projection, and the subscript i means that the operator acts on 
the i'th component in the A^-fold tensor product appropriate for A^-particle 
systems. Working with the right side of (0) is often easier than handling the 
left side, using the fact that it is a superposition of products of one-particle 
projection operators. 

In Section ^ below we will use the decomposition of the Yukawa poten- 
tial to derive the high density asymptotics of the ground state energy of 
jellium with Yukawa interaction in the thermodynamic limit. This is a gen- 
eralization of [|GS| , where an analogous expansion in the case of Coulomb 
interaction was accomplished. 

Yukawa potentials are actually used in solid state physics as a model for 
the screened Coulomb interaction between electrons in a solid (see, e.g.. 



2. The Decomposition 

In this section, we state and prove the decomposition (y), and comment 
on its implications. In the following theorem, we restrict ourselves to the 
case n > 2, and remark on the easy case n = 1 after the proof. 

THEOREM 1 (Decomposition of Radial Functions). For n > 2, let 

V : M" — > M 6e a radial function that is [n/2] + 2 times dijjerentiable away 
from X = 0. For m G Nq denote V^"^\\x\) = d"^ / d\x\"^V (x) . Assume that 
lim|^,|^oo kri^^™H|x|) = for all < m < [n/2] + 1, and let Xr{x) = 
9{r — \x\). Then 



V{x)= drg{r)xr/2*Xr/2ix) , (3) 
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where 

(_l)[n/2] 2 



g[r) 



T{'^) (7rr2)("-i)/2 

/ dsy(["/2l+2)(,)M j ,(,2_,2)2(n-3) 

+ <5odd^""/'l+')(r)r(2r)^(-3)r(2i^)) , (4) 



and 5odd = 1 /c"^ n odd, 5odd = for n even. 

Note the r/2 in (|3|), which is chosen for convenience. 

Proof. Elementary considerations show that 



1 /7r\("-l)/2 r 2 2^i( 

Xr/2*Xr/2(3^) = ^^^^^ (^^J J dy{r - y )2^ 



(5) 



for l^l < r, and otherwise. Inserting the definition (^) for (7, we can 
therefore write 

l)[n/2] 1 ^~\ 

drg{r)xr/2*Xr/2{x) 



~,n—l I ^ ' 



(J \ 1— 1— [n/2] 1 /■'■ 1 

^ j S{S^ - r2)2("~3) / ^y(^2 _ y2)2(n-l) 

dr_L_y(K2]+2)(^)^(2r)2("-3)r(i^) 
X r dy{r^ - y^)^^''-^^ . (6) 



We now use the fact that 

r 1/2 2\o(n-l) / 2 2\TT(n-3) 

dr — 3T(r — y j2^ 's[s —r j2^ ^ 



^ ^ n-l-[n/2] 



ds / Jy r" ^ 



n-l-[n/2] 



+ i<.dd ^(." - 9=)l'"-'>.(2s)5l"-'lr(Sjl) , (7) 
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and change the order of integration to get 

(_l)["/2] I \^ roo 

drg{r)xr/2*Xr/2{x) 



r(n_i)r(H±i)2-2y ^0 

dsy(["/2]+2)(^) / dy (8) 

x\ ^Ix 

, \ n-l-ln/2] />s 1 1 1 

X ( - ) / dr^s{s' - r2)2("-3)(r2 - y2)2("-i) . 

dsj Jy r" ^ 



The last integral can be evaluated to be 



(9) 



sT{n) 
and therefore 

drg{r)xr/2 * Xr/2{x) 

r([n/2] + 2) 4| ^ ^^ ' '^ V ;, V ; 

where we integrated by parts in the last step, using the demanded decrease 
properties of |xpy*^™)(|x|). D 



f) holds with g{r) = V (r). This is a well known result 
, and was used in ||HS[ to get an estimate on the indirect 



For n = 1, ( 
due to Polya 
interaction energy in one-dimensional quantum system. 

REMARK 1. If V{x) is n + 1 times differentiable away from x = 0, we 
can use partial integration to write g{r) in the more compact form 

/'_i\n+l 9 foo 

y^ > r(^) (7rr2)("-i)/2 J^ "■ ' "■ ' "■ ' 

For simplicity, we have restricted ourselves to formulating Theorem [llfor 
differentiable functions only, but as is obvious from the proof, (^) and (11) 
hold for a more general class of potentials, if the derivatives are interpreted 
in the sense of distributions. 

EXAMPLE 1. A particular simple example which our decomposition ap- 
plies to is the Coulomb potential V{x) = l/|x|. For the general case of n 
dimensions we compute the weight function to be 

, , I Tin + 2) 7ri-"/2 
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EXAMPLE 2. Another example, which will be used below, is the Yukawa 
potential Yfj_{x) = exp{—fi\x\)/\x\, with /i > 0. In three dimensions, the 
corresponding weight function reads 

gir) = -^:^ [8 + 8fir + i{f,rf + ifirf] . (13) 

The decomposition @ provides conditions for positive definiteness of V. 
By positive definiteness of a locally integrable function V we mean that 



V{x — y)ip{x)ip{y)dxdy > (14) 

for all ip G Co°(M"). This is equivalent to V having a positive Fourier 
transform in the distributional sense. Since Xr * Xr{x) is obviously positive 
definite (as a function of x), any V with a positive weight function g{r) is 
positive definite. It is easy to see that 

(r \ n— 1 — [?i/2] -j^ 

— J s(s2-r2)2("-3) >o for s>r, (15) 

and therefore positivity of (— l)["/2]y(["/2]+2) jniphes that V{x) is positive 



definite. This is a result due to Askey [Al], who deduces it from certain 
positivity properties of Bessel functions proven in |A2] and |p^ . Further 
extensions and refinements can be found in [ p2| and references therein. 

Consider now the (physically interesting) special case n = 3. Here the 
weight function g takes the simple form 

air) = —, (F"(r) - rV"'{r)) = -^ (^^] . (16) 

Hence a monotone decrease of V /r (instead of y ) is sufficient for positive 
definiteness of V. This condition for positive definiteness of a radial function 
on M.^ has been obtained before by A. Martin [GM|. See also [HN] for an 
earlier reference. 

Positivity of g might be a useful way of checking positive definiteness of 
a function V. However, this condition is of course not necessary, as the 
following example shows. 

EXAMPLE 3. For n = 3, V{x) = l/(|xp + 1) is positive definite, but the 
corresponding g{r) is not a positive function. 

3. Ground State Energy of Yukawa Jellium 

We now demonstrate the usefulness of our decomposition (|3|) in applica- 
tions to A'-body quantum systems. We consider a model of electrons in- 
teracting with Yukawa potentials and with a uniform charged background. 
Yukawa potentials are used in solid state physics as a model for screened 
Coulomb potentials Q. They were first introduced by Debye and Hiickel 



[DH], and were later used in meson theory by Yukawa. 
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The system under consideration is described by the Hamiltonian 

N 
j=l l<i<j<N 

acting on 7^ = Aj=i ^'^{^^ dxi] C^), the antisymmetric tensor product of the 
one-particle space L^(A, dx,; O) of one electron with spin {q — l)/2. The 
particles are confined to the cube A = [0, L]^ of side length L, with volume 
|A| = L^. They interact via the Yukawa potential 

p-fj.\x\ 

Y,{x) = —- , (18) 

\x\ 

with yU > 0. The electrons move in the potential V created by a charged 
background, given by 



V{x) = p dyY,,{x-y) , (19) 

Ja 

where p > is the background density. We assume neutrality of the system, 

i.e., p = iV/|A|. 

We have to specify boundary conditions for H on A, and we take Dirichlet 

boundary conditions for simplicity. Denoting E{N, A, /x) = infspec H, we 

define the energy density for a neutral system in the thermodynamic limit 

as 



lim — ( E(p[A[, A,/i) + ip2 / dxdyY^ 

L^oo |A| V JAxA 



e{p,^l) = lim — E{p\K\,K,p) + \p^ / dxdyY^{x -y)] . (20) 



For convenience, we have added the self energy of the background charge. 
For p = this is essential, because otherwise the limit does not exist. We 
take the existence of the limit (|20| ) for granted. For the Coulomb case (i.e. 
^ = 0), existence (and independence of the boundary conditions for H) of the 
thermodynamic limit has been shown in [|LN| . However, because of the short- 
range nature of the Yukawa potential, the existence of the thermodynamic 
limit for p > is much easier to show. 

Our main theorem concerns the behavior of e(p, p) for large p. This is 
an extension of the result in [G£] for p = 0. We denote by kp the Fermi 



momentum 

kF = {Gn^p/qy/' . (21) 

THEOREM 2 (Ground State Energy of Yukawa Jellium). As p ^ 

CX), 



e{p,p) 



uniformly in p for bounded p/kp. Here J is the function 

J{v) = ^-l'i^ + lyi^ + 1^^') 1" (i + ^) - ^^ ^'^cta'^ I ■ (23) 
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As in [ pS(| , 0(1) is positive and smaller than 0{p^^'^^^'^) for any e > 
and for fi/kp fixed. The second term on the right side of ( [2^ ) is entirely due 
to correlations in the ground state wave function and is usually referred to 
as "exchange energy". The function J{rj) in Theorem |^ is in fact the integral 



Jiv) 



dk 



e 



k^ + r/^ 



1 



-k-\ k-^ 

I 16 



(24) 



Notice that J(0) = 1, whereas J{r]) « 4/(97/^) as rj —>■ 00. 

We conjecture that the asymptotics in ( p2|) is uniform in /i, even if fi/kp 
tends to infinity. For q = 1, i.e., the spinless case, this is in fact easy to 
prove, as will be remarked after the proof of Theorem |2|. Note, however, 
that in physical situation n/kp is small. The range of the potential is given 
by 1/fi, which is always bigger than the mean particle distance p^^'"^. 

The proof of Theorem |2| follows essentially the analogous discussion of 
the Coulomb case in jGS| . We use freely the estimates derived there, mainly 
pointing out the differences to our case. 

Proof. To obtain an upper bound to the ground state energy, we use as a 
trial state a Slater determinant of the A^ lowest eigenvectors of the Laplacian 
on A, including spin. The calculation of the expectation value is essentially 
the same as the one done by Dirac [p| in the Coulomb case. In the thermo- 



dynamic limit the boundary conditions do not matter (see [|GS[1 for details), 
and one can just consider plain waves, with momenta up to the Fermi mo- 
mentum kp. The function J then arises from the integral 



dk 



(2vr) 



y,ik) 



3"^ 



dq 

R3(2^ 



eikp - \q\)eikF - \k - q\) 



i^^^H^ 



(25) 



We are left with the lower bound. We start with the following decom- 
position of the Yukawa potential, proven in Section ^. If Xr denotes the 
characteristic function of a ball of radius r centered at the origin, then 



where 



y,i{x-y) 



g{r) = 2g{2r) 



dz I drg{r)xr{x - z)xr{y - z) , 
/o 



1 6-2/^*^ 



vr 



[l + 2pr + 2{prf + {prf 



(26) 



(27) 



(compare with Exam ple || ). The key to the lower bound is a correlation 
inequality derived in |GS, Cor. 5]. It states that for any projection X and 
< P < 1 on the one-particle space, and for antisymmetric ^ € 7Y, 



^ 



zZ -^^-^i 



l<i<j<N 



^ ) > iTr[A7]2 - ^Tt[PXPX] 



-const.Tr[A(P + 7)]min|l,(Tr[A(l-P)7(l-P)])^/n . (28) 
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Here we denote by 7 the one-particle reduced density matrix of ^|J, with 
corresponding density p^. Using this, with X = Xri' ~ z), and integrating 
over r and z as in (^), we get that 



i; 



7 ^ JL^[Xi Xj^ 
l<i<j<N 



^) >\ j dxdyp^{x)p^{y)Yfj_{x - y) 
/ JAxA 



'\^ I dxdy\P{x,a]y,(T')\^Y^{x-y)-^Tioi . 
I J AxA 



(29) 



The function P{x,a;y,a') is the integral kernel of P, and a € (l,...,g) 
denotes the spin variables of the electrons. The error term is 



Error 



const. dz drg{r)[xr*{P'r + Pp)]iz)inm\l,[xr* PQ^QJiz)^^"^} ■ 

JR3 Jo ^ -" 



(30) 



Here we introduced the operator Q = 1 — P; pp and pq^q are the densities 
corresponding to P and Q^Q, respectively. 

Since g{r) < const. r~^, we can proceed exactly as in |GS| , Lem. 6] to 
estimate 



Error < const. II p^ +PP 11]^' ^IIPt + PPII5/3 '^(7) -P)"*" "^ 



(31) 



for any < e < 1/6. Here 5{'y,P) = Tr[7(l — P)] measures the "difference" 
of 7 and P. Since ip is supposed to be the ground state, we can use the 
upper bound to get the a priori knowledge that II/O7IL/3 < const. A^/cp, see 



5/3 



m, Eq. (4.11)]. 

Now let P be the projection onto the first A^ eigenstates of the Laplacian 



with periodic boundary conditions. It is shown in |GS, Eq. (4.13)] that, for 
any Tp G H, 



N \^ 



N 






A>4(^ + const. (iV-^<^(7, P)) ^] - 0(1 



) (32) 



as |A| -^ 00. The second term in (P9|), when divided by N, converges in 
the thermodynamic limit to ■^kpj{p/kp), as explained in the upper bound. 

1 5/3 



Moreover, II/OPII5/3 < const. Nkp. Putting together (p9|), (p|) and (^), we 
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therefore get that in the ground state ip, 

^ (V' 1^1 V') + ^P' y dxdyY^ix-y) 

> 4 ("^ + const. {N~'6{j,P)f^ - ^kFJifi/kp) 

- const. A;f {N-^6{j, P)) ^^^'' 
+ \ I dxdy{p^{x)- p){p^{y)- p)Y^{x-y)-o{l) (33) 



'AxA 

as |A| — > oo. The last term is positive, since 1^ is positive definite. Mini- 
mizing the right side of (^) over (5(7, P) gives the desired result. D 



REMARK 2. We expect the asymptotics in (22) to be uniform in p/k^ 
even if this value goes to infinity. The reason is that for p/k^ ^ 1 the in- 
teraction potential is of very short range, and hence the interaction between 
particles of the same spin can be neglected to leading order due to the anti- 
symmetry of the wave function. I.e., in this case one expects a contribution 
to e(/9, p) from the interaction energy of (1 — q~^)2'K p^ / p^ . The total energy 
for p » p^i^ should then be 

e(/3,/x) 3 2 47r , _^ 27r 27r 

which is exactly the same as (22) as p/k-p -^ oo. 



For the spinless case, i.e., q = 1, the uniformity is in fact easy to prove. 
Neglecting the positive interaction energy we get a lower bound 

e(Pi/^) ^32 27r 

> 7% J/' ' 35 

p b p'^ 

which agrees with (p2|) as p/k^ — > oo. We note that the second term is now 
not the exchange energy, but the difference of the energy of the electrons in 
the background and the self energy of the background. 
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